Universality and robustness of revivals in the transverse field XY model 



Juho Happola, 1 " 2 ' 3 Gabor B. Halasz, 3 ' 4 and Alioscia Hamma 3 

1 Aalto University, Helsinki, Finland 
2 University of Helsinki, Finland 
3 Perimeter Institute for Theoretical Physics, 31 Caroline St. N, N2L 2Y5, Waterloo ON, Canada 
4 Trinity College, University of Cambridge, Trinity Street, Cambridge CB2 ITQ, UK 

We study the structure of the revivals in an integrable quantum many-body system, the transverse field XY 
spin chain, after a quantum quench. The time evolutions of the Loschmidt echo, the magnetization, and the 
single spin entanglement entropy are calculated. We find that the revival times for all of these observables are 
given by integer multiples of T rcv ~ L/u max where L is the linear size of the system and t> m ax is the maximal 
group velocity of quasiparticles. This revival structure is universal in the sense that it does not depend on the 
initial state and the size of the quench. Applying non-integrable perturbations to the XY model, we observe 
that the revivals are robust against such perturbations: they are still visible at time scales much larger than the 
quasiparticle lifetime. We therefore propose a generic connection between the revival structure and the locality 
of the dynamics, where the quasiparticle speed u ma x generalizes into the Lieb— Robinson speed vlr. 



PACS numbers: 

I. INTRODUCTION 

The behavior of quantum many-body systems away from 
equilibrium has recently become the object of intense experi- 
mental study. Ultra-cold atoms in optical lattices feature both 
large phase coherence times and a high degree of controllabil- 
ity, making it possible to observe quantum coherent dynamics 
1 1 1. The potential technological implications are profound be- 
cause the manipulation of coherent quantum dynamics is at 
the root of the possibility of building a quantum computer. 

In terms of theoretical description, systems away from equi- 
librium are considerably more complicated than their equi- 
librium counterparts. Whereas equilibrium systems can be 
understood by means of standard methods like mean field 
theory and renormalization group, we lack analogous meth- 
ods for understanding non-equilibrium physics \2\. There are 
no obvious generalizations of such standard methods at equi- 
librium, and in particular, it is not clear to what extent the 
non-equilibrium dynamics of quantum systems features uni- 
versality. An important example is the universality of the 
Kibble— Zurek mechanism to compute the density of defects 
as the external temperature is tuned in time (3) . 

Experiments on cold atomic gases give us the opportunity 
to observe a genuine quantum evolution in a system that is 
very close to be isolated. We can then address questions like 
the possibility of equilibration in a system with unitary evo- 
lution and the process of thermalization in a closed quantum 
system. These topics have recently found a renewed interest 
along with other fundamental problems in quantum statisti- 
cal mechanics (U[5)- In the case of such systems, the non- 
equilibrium dynamics is obtained by making some parameters 
of the system Hamiltonian time-dependent and hence cast- 
ing the system off equilibrium. The time dependence of the 
Hamiltonian can be adiabatically slow [6| or it can change 
abruptly Q O in which case the process is called a quantum 
quench. Among other applications, the paradigm of the quan- 
tum quench has been recently used to study the behavior of 
topological phases away from equilibrium (9j. 

The understanding of non-equilibrium dynamics of isolated 



interacting quantum systems is of fundamental importance to 
understand quantum equilibration, a topic that has lately ex- 
perienced a new renaissance (4j [5] [TO] . Due to unitary time 
evolution, a finite quantum system with a non-trivial initial 
state can not converge to a steady state. However, it has been 
shown in some remarkable papers lfTTl[T2ll that thermalization 
of a finite subsystem is possible in the infinite size limit. The 
system can be coarse grained by choosing a partial set of lo- 
cal or macroscopic observables, and the expectation values of 
these observables can in principle converge to the ones in ther- 
mal equilibrium. In particular, an integrable system does not 
thermalize: even local or macroscopic observables undergo 
oscillations, or at best relax to equilibrium values that are gen- 
erally not the same as those predicted by the microcanonical 
ensemble lfT3l FBI . Even for a non-integrable system, ther- 
malization does not always occur: in some cases there is only 
relaxation to a non-thermal state depending on the initial con- 
ditions irrsi . 

The equilibration process has several characteristic time 
scales. The largest one is the recurrence time T rcc at which 
the system gets infinitesimally close to its initial state. It only 
exists for finite systems, and diverges at least exponentially as 
a function of the system size. The smallest time scale is the 
relaxation time T Te \ at which a given observable relaxes to its 
long-term average value. Furthermore, there is a third impor- 
tant time scale T rov in between at which revivals occur: these 
are brief detachments from the average value of an observable. 
Revivals typically last for a very short time only (in compari- 
son to the spacing between them), and their magnitude decays 
in time as the equilibration process nears completion. 

In this paper, we investigate the structure of the revivals for 
an integrable model, the transverse field XY spin chain. Our 
main result is that this structure is universal in the sense that 
it does not depend on the details of the quench and on the ini- 
tial state. By applying non-integrable perturbations to the sys- 
tem and finding that the revival structure is surprisingly robust 
against such perturbations, we argue that the revival structure 
is a universal non-equilibrium property which follows from 
the locality of the Hamiltonian. 
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II. FORMALISM OF THE QUANTUM QUENCH 

We consider a closed quantum system whose Hamiltonian 
if (Ai, .., Xr) depends on parameters Xi representing the cou- 
pling strengths of interactions and external fields. A quan- 
tum quench is a sudden change in the Hamiltonian of the 
system. The quantum system is originally prepared in the 
ground state p of H(X^), and at time t = we switch 
the parameters to different values X^ 2 \ The system then 
evolves unitarily with the quench Hamiltonian H(X^) ac- 
cording to p(t) = Ut{po), where we define the superoperator 
U t (X) = exp(— iHt)X exp(iHt). Unitary evolution implies 
that a finite system can not converge to a steady state p, even 
weakly. The limit of p(t) for t —> oo does not exist unless 
the initial state po is trivial, e.g., an eigenstate. On the other 
hand, the time average p = limbec t~ 1 J Q p(s)ds always ex- 
ists and is given by the po totally dephased in the eigenbasis 
n„ = \E n ) (E n \ of the Hamiltonian: p = J2 n n n /J n n . For 
a finite system, equilibration means that the expectation val- 
ues of macroscopic observables spend most of their time very 
close to their average values. 

An important quantity describing the time evolution is the 
Loschmidt echo (LE) defined as C(t) = |tr(exp(— itH)po)\ 2 
which gives a measure of the distance between the time 
evolved state p(t) and the initial state p . When the system 
undergoes a recurrence at t — T ICC we have C(T TCC ) ~ 1. The 
general expression for the LE can be written as 

C(t) = J2PnPme- i ^- E ^ t , (1) 

n,m 

where p n are the populations of the Hamiltonian eigenstates 
for the initial state. It follows that the time average of the 
LE is C = tr(p 2 ) = ^2 n Pn- The LE typically decays in a 
short time T TC \ from 1 to its average value C around which it 
oscillates. The relaxation time T re i is O(l) for an off-critical 
quench, while it scales like O(L^) with the system size L for 
a critical quench, that is, if A^ 2 ' is a critical point, showing a 
critical slow-down of the system lfl6l[T7l . 

Revivals are also visible in the LE as deviations from the av- 
erage value C. We define revivals as time instances at which 
the signal C(t) differs from L by more than three standard 
deviations. According to Eq. Q, this happens when an ex- 
ceptionally large number of weights p n get partially back in 
phase. It is not straightforward to understand directly from 
Eq. ([T]i when such a situation can occur in a generic quantum 
system, therefore we consider the particular case of a simple 
one dimensional spin chain. 

III. REVIVALS IN THE XY MODEL 

A. Exact solution by free fermions 

In this section, we consider an integrable (exactly solvable) 
model, the one dimensional XY model of N spins one half in 
a transverse magnetic field. Since the model is exactly solv- 
able, we can obtain the whole spectrum and the eigenstate 



decomposition of the initial state. This leads to an exact ex- 
pression for the LE, and the revival times can be extracted by 
inspecting its time dependence. 

The Hamiltonian of this spin chain is given by 

N , 

i=i ^ ' ' 

(2) 

where 77 is the anisotropy parameter, and h is the external 
transverse magnetic field. We assume cyclic boundary con- 
ditions (Tn+i = (— l) 9 cri. The periodic (q = 0) and antiperi- 
odic (q = 1) boundary conditions differ in 0(l/N) terms and 
this difference usually does not affect the phase diagram or 
other quantities in the thermodynamic limit. However, it can 
be important in the LE that is typically exponentially small in 
TV. In the following, we shall see that the boundary conditions 
can have a dramatic effect in the case of the critical quench. 
Note that the XY model reduces to the quantum Ising model 
for 77 = 1, and to the isotropic XX model for 77 = 0. The 
Hamiltonian exhibits two regions of criticality: the XX model 
at 77 = has a critical region for h E (—1,1), while the XY 
regions of criticality are the lines h = ±1. 

The XY model can be diagonalized by a standard procedure 
lfT8ll . In the first step, we map ai to spinless fermions by using 
a Jordan— Wigner transformation: 

of = 1 - 2c\c h of = ( ( r+) t = cfe i7r ^=i c h . (3) 

The translational symmetry is then exploited by the Fourier 
transform c; = YI1L1 elklc k, where the momenta are 
quantized according to k„ = tt(2ti + 1 — q)/N. Finally, after 
the Bogoliubov transformation Cf. = cos 6u1k + i^^^kl^^ 
we obtain a Hamiltonian describing non-interacting fermionic 
degrees of freedom jk- 

k>0 

The dispersion relation of these fermionic quasiparticles is 

given by Afc = ^/e 2 + rf- sin 2 k with ek = h — cos k, and 
the angle Of. appearing in the Bogoliubov transformation is 

Ok = tan- 1 [(?7sin£;)/(e fc + A fc )]. 

B. Loschmidt echo and revival times 

In a quantum quench, the system is prepared in the ground 
state of 7J(A (1) ), and then evolved with ff(A (2) ) at t > 0. It 
is useful to write the ground state \ip (0)) of the initial Hamil- 
tonian in terms of the eigenstates of the quench Hamiltonian: 

l^(0)> = Y[(cosxk-ism X klW-k) I *) » ( 5 > 

k>0 

where \k = 0^ ~ ^fc^' an< ^ l^fc) is the vacuum state defined 
by jk |0fe) = 7_fc |0fe) = 0. The time evolution then reads 
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\i>(t)) = exp(-iHt) \ip (0)), and the LE takes the form 

c w = n ( i - Aks[n2 ( Afc *)) - 

fc>0 



(6) 



where the coefficient Ak = sin 2 {2\k) is a slowly varying 
positive function of the momentum k. 

Now we show how the revivals in the LE can be derived 
from the dispersion relation of the quasiparticles. We first 
take the logarithm of Eq. |6} to transform the product into a 
sum over momentum k: 



ln£(i) = ^ln(l-^ fc sin 2 (A fc t)) 



(7) 



k>0 



Since each k > term has a periodicity ir/Ak in time, and 
Ak varies sufficiently slowly with k, this expression shows 
that nearby k modes separated by A/c = 2ir/N add up con- 
structively whenever AA^i = pir with p E Z. This rear- 
ranges to tk = \pN \dKk/dk\ , and in principle we could 
expect a revival time tf. corresponding to each k. However, 
more modes can add up constructively if the dispersion re- 
lation Aft is closer to a straight line, therefore the most pro- 
nounced revivals (in fact, the only revivals that stand out from 
the background noise) are given by the stationary values of 
the group velocity v g (k) = \dh.k/dk\. The first of these re- 
vivals is the one corresponding to the maximal group velocity 
u mM = maxfe v g {k) and we can thus give the following esti- 
mate for the revival time scale: 



N 



2v n 



N 
~2 



Ok 



(8) 



The maximum group velocity w max can be computed exactly 
from the dispersion relation A*.. It takes a particularly simple 
form in the case of the Ising model (j) = 1): v max = h when 
h < 1 and v max = 1 when h > 1. 

In Fig. 1(a) the LE for a critical quench (/12 = 1) with 



antiperiodic boundary conditions (q = 1) is plotted. The T rev 
predicted by Eq. (|8]l is in perfect agreement with the data. 
The same critical quench with periodic boundary conditions 
(q = 0) is interesting. At the odd revivals, there is no signal 
in the LE due to a destructive (vanishing) contribution from 
one of the k modes in the product of Eq. (|6j. Only the even 
revivals are spotted. 

Figs. 1(b) and 1(d) verify that the revival times scale like 
0(L) where L ~ N is the linear size of the system. The 



quench in Fig. 1(b) is non-critical, and the parameters of the 



quench Hamiltonian are far away from any phase boundaries. 
This means that the system is gapped, and a simple spectral 
analysis would imply that T rcv is of O(l). However, even if 
there is a small reconstruction of the wave-function at a time 
scale 1/ A (where A is a difference between any two energy 
levels), the weight involved is not sufficiently large to make 
the corresponding revival strong enough. Visible revivals are 
governed by Eq. dHJ. 
In Fig. 

As 



1(c) the system is quenched from different ground 
states corresponding to different parameter values 




one can see, the details of the evolution and the average values 
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FIG. 1: Loschmidt echo for the quenched XY model with an- 
tiperiodic boundary conditions (q = 1). |(a)| Critical quench with 
Vi =T] 2 = 2.0, hi = 0.5, h 2 = 1.0, N = 400, and u max = 2.0. 
Eq. l[8j gives T rcv = 100.0, and from the plotted data we get 
T rC v ~ 99.7. I (b)| Non-critical quenches for different system sizes: 
Vi = m = 2.0, hi = 0.7, hi = 0.8, N € {200,400,800, 1600}, 
and iWx ~ 1.77. The equilibrium value of the LE decreases with 
system size. Eq. |8| gives T rcv ~ 112.7 for iV = 400, and from 
the data we get T rev w 113.0. |(c)| Quenches of different size: 
Vi = 772 = 2.0, hi G {0.5,0.6,0.7}, h 2 = 0.8, and N = 400. 
Group velocities and revival times match those in subfigure |(b)| |(d)| 
Linear scaling of the revival times with system size. The lines rep- 
resent the values given by Eq. {8}. Parameters used: 1: those in 
subfigure |(a)| 2: those in subfigure |(b)| In each plot, the red hori- 
zontal lines represent average values and three standard deviations 
thereof. 



C are different, but the structure of the revivals is the same 
for all the quenches, and the revival times are consistent with 
those predicted by Eq. (JHJ. This is the promised universality 
of the revival structure: the initial state and the size of the 
quench are unimportant. The parameters A 1 - 2 ' determine i> max 
and therefore T rov but not the fact that revivals happen at time 
instances spaced evenly at intervals that are linear in system 
size: t p = pL/v max . 



C. Magnetization and entanglement entropy 

Although the behavior of the LE is illuminating for theoret- 
ical arguments, it is hardly of experimental relevance because 
the amplitude of the signal is zero in the thermodynamic limit. 
On the other hand, it is expected that a revival in the LE cor- 
responds to revivals in macroscopic observables as well. To 
demonstrate this, we study the time evolution of the order pa- 
rameter (magnetization) fi(t) ETl and of the single spin entan- 
glement entropy S(t). If we consider a subsystem A consist- 
ing of a single spin, the von Neumann entropy between sub- 
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system A and the rest of the system is S — — Tr (pa log pa), 
where pa is the reduced density matrix of subsystem A. 
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FIG. 2: A small quench with q = 1, 771 = 7)2 = 2.0, /ii = 0.7, /12 = 
0.8, AT = 200, and u ma x ~ 1-77. From top to bottom: Loschmidt 
echo C(t), magnetization n(t), and single spin entanglement entropy 
S(t). Eq. |8} gives T rev ~ 56.37, and from the plotted data we get 
T rov ~ 56.58. The horizontal lines represent average values and 
three standard deviations thereof. 

When the initial state of the quenching process in Eq. <[5j is 
expanded in position basis rather than momentum basis, one 
finds that the excitations are pairwise correlated between dif- 
ferent lattice sites due to the fermionic anticommutation rela- 
tions. This implies that pa is diagonal, and then we can use 
translational invariance to establish 

5 (t) = -11 (t) \ogp (t) (t)) log (1-p (<)), (9) 

where the average magnetization is given by 

1 b 



(10) 



{c\c k ) = sin 2 (6» fe ) cos 2 (xk) + cos 2 (6 k ) sin 2 (xfe) 

- 2 sin (6k) cos (9 k ) sin (xk) cos (xk) cos (2A k t) . 

The entanglement entropy S and the average magnetization p 
are thus governed by the interference of the same modes as in 
Eq. |6} with their frequencies given by the same dispersion 
relation A*.. This means that they share the same time scales 
as the LE. The results are illustrated in Fig. [2] In fact, since 
the magnetization is given as a sum rather than a product of 
different oscillating modes it can be used instead of the LE to 
determine T rev in the large N limit when the LE has a very 
small average value. In practice, the magnetization is analo- 
gous to the Logarithmic Loschmidt Echo used in lfT9l . 



D. Time evolution of a local disturbance 

Now we investigate the time evolution after a local distur- 
bance in the spin chain. In particular, we consider a single 
spin flip at position I, which is represented by the operator 
F[ = c\ + cj. The time evolution of the resulting local dis- 
turbance is best studied in the Heisenberg picture, where the 



operator Fi becomes time dependent and takes the form 

Fi(t) = [ni-i>(t)ci> + nu,(t)4] , (li) 



fil-J'(*) 



I' 

1 

N 



,ik(l-V) 



e 4Afc *sin 2 9 k + e- lAkt cos 2 6 k 



+ i(e 



sin u k cos Ok 



This expression shows that nearby modes at momentum k add 
up constructively whenever AA k t — 2pn ± Ak\l — l'\ with 
p € Z. Due to Ak = 2ir/N, this condition can be written as 
t= (pN ±\l- l'\) \dA k /dk\~\ and one verifies that the dis- 
turbance travels with the group velocity v g (k) — \dA k /dk\. 
Once again, more such modes can add up constructively if the 
second derivative of A k vanishes, therefore we expect visi- 
ble wave packets to travel with the stationary values of v g (k). 



Indeed, Fig. 3(a) shows that the different wave packets corre- 



sponding to the local extrema of v g (k) propagate through the 
lattice while maintaining their respective wave forms. 
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FIG. 3: |(a)| Time evolution of a local disturbance (local quench) 
\Qi(t)\ 2 with q = 1, tii = 772 = 2.0, hi = 0.7, h 2 = 0.8, and 
N = 200. The maximum speed is ?; max w 1.77. The corresponding 
global quench is illustrated in Fig. |l(b)| for comparison. |(b)| Maximal 
group velocity fj max as a function of 77 and h. 

The fastest wave packets travel with the maximal group ve- 
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locity t> max , and the first revivals can be interpreted as con- 
structive interferences between them. Since we assume cyclic 
boundary conditions, one can think of the spin chain as a 
closed ring. In this picture, the fastest wave packets first meet 
halfway in the ring at time t = N/2v max which indeed coin- 
cides with the first revival time in Eq. (IS). 



IV. NON-INTEGRABLE PERTURBATIONS 

In the previous section, we showed that the structure of the 
revivals is governed by the maximal speed of quasiparticles in 
the system. We found that as long as the quasiparticles exist, 
the details of the quench are not relevant. There is universal- 
ity within the integrable behavior of the system. At this point, 
we wonder whether there is universality beyond the quasi-free 
system. The local physics induced by the local Hamiltonian 
might imply that as long as information is not completely lost 
— as in the case of an infinite system — revivals can be ob- 
served due to the recombination of the fastest signals even 
when these are not point-like and do not correspond to quasi- 
particles. To investigate this possibility, we now study the ro- 
bustness of the revival structure against non-integrable pertur- 
bations. 



A. XZ spin chain 

We start by considering the XZ spin chain, which contains 
an additional of of +1 coupling with respect to the quantum 
Ising model (77 = 1): 



H 



1 N 



f of+i + 90' of+i + hof 



(12) 



1=1 



This Hamiltonian is non-integrable, and we simulate the quan- 
tum quench by exact numerical diagonalization. To achieve a 
relatively large system size (N ~ 50), we take the limit of 
large field (h ^> 1), restricting the effective Hilbert space to 
states where almost all spins are aligned with the field. 





(a) 



(b) 



FIG. 4: Loschmidt echo for the quenched XZ spin chain with peri- 
odic boundary conditions in the integrable case g% = (a) and the 
non-integrable case g% = 0.3 (b). The other quench parameters are 
r?i = 772 = 1.0, hi = 10.0, h 2 = 4.0, and gi = 0.0 in both sub- 
figures. Since N = 40 and t> ma x = 1, Eq. ([HI gives T rcv = 20. 



As shown in Fig. [4] the structure of the revivals is clearly 
visible for \g\ < 0.5, and hence this structure is universal for a 
range of the non-integrability parameter g. The revival times 
scale linearly with N, and the range of visibility is largely in- 
dependent of both the system size N and the magnetic field h. 
We find that pronounced revivals gradually disappear in the 
range 0.3 < \g\ < 0.7 for all 20 < N < 50 when h = 4 
is fixed, and for all 2 < h < 100 when N = 20 is fixed. 
Since v max ~ 1 in the Ising model for h > 1, we see that the 
non-equilibrium dynamics is dominated by the afaf +l term 



in Eq. ( 12 1. The XZ spin chain is therefore significantly non- 
integrable for \g\ ~ 0.5, even in the h — > 00 limit. This claim 
is further supported by the fact that the visibility range in g 
does not depend on h. We finally note that the equilibrium 
(long-term average) value of the magnetization is strongly de- 
pendent on the initial state and so this equilibration is not ther- 
malization, even though the system is non-integrable 1 15 1. 



B. Random disorder in the field 

Now we consider another integrability breaking perturba- 
tion to the XY spin chain. We introduce a site dependent ex- 
ternal field component e; to the Hamiltonian that explicitly 
breaks the translational invariance and the integrability of the 
model. The total site dependent field is hi = h + e/,, where ei 
is randomly picked from a uniform distribution with a maxi- 
mum amplitude e in the sense — e < ej < e. 

Using exact diagonalization and exploiting the fact that the 
Hamiltonian decomposes to subspaces of odd and even num- 
ber of spins, we can assess the effect of the site dependent 
field disturbance for modest sized systems N < 13. The sim- 
ulations run on such systems give supporting evidence that 
the revival structure is essentially unchanged for a range of 
the amplitude e. In Fig. [5] the Loschmidt echo for a critical 
and a non-critical system is plotted with various values of e to 
illustrate this observation. 



V. DISCUSSION 

We found in SectionUUthat the transverse field XY model 
exhibits a universal structure of revivals that is independent of 
the initial state and the size of the quench. Since the XY model 
is integrable and its exact solution is in terms of free fermionic 
quasiparticles, there is a straightforward interpretation for the 
phenomenon of revivals. The information propagates around 
the system via wave packets of quasiparticles, and the first 
revival occurs when the wave packets traveling with the maxi- 
mal group velocity w max meet. This interpretation explains the 
universality of the revival structure since i> ma x only depends 
on the dispersion relation of the quasiparticles associated with 
the quench Hamiltonian (and nothing which would be related 
to the initial state). 

On the other hand, the robustness against non-integrable 



perturbations found in Section IV suggests something more 
generic than the quasiparticle interpretation. The quasiparti- 
cles of the perturbed system are no longer free but there is a 
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FIG. 5: Loschmidt echo for the transverse field Ising model when 
the external field has a small site dependent component of varying 
magnitude: 771 = 772 = 1.0, N — 13, periodic boundary conditions. 

f Non-critical case: hi = 0.5, hi = 0.8, ei = 0.0, £2 = e. Eq. 
gives T rcv = 8.12 for e = 0.0. [(b)] Critical case: hi = 0.5, 
h 2 = 1.0, ei = 0.0, e 2 = e. Eq. l|8j gives T rov = 6.5 for e = 0.0. 

finite interaction between them. For the XZ spin chain, this 
interaction is on the order of g as can be verified by look- 
ing at the exact eigenvalues for finite systems. This implies 
that the quasiparticles decay on the time scale of l/g, and 
hence one would not expect to see revivals at time scales much 
larger than l/g. However, the evidence is on the contrary: 
at N = 50 and g ~ 0.5, revivals are still clearly visible at 
t ~ 50, which is an order of magnitude larger than l/g ~ 2. 

It appears that the revival structure is a non-equilibrium 
property that is beyond integrability and the existence of sta- 
ble quasiparticles. Here we provide a more generic interpre- 
tation in terms of locality. Quantum many-body physics is 
generally described by Hamiltonians that can be written as 
sums of local operators. The locality of the Hamiltonian has 
profound consequences on the dynamics of the system ll20l : 
there exists an emergent light cone for the propagation of in- 
formation such that signals outside the light cone are exponen- 
tially suppressed. The characteristic speed of the light cone 
gives the maximal speed of information in the system, and it 



is called the Lieb— Robinson speed v^r. In general, it de- 
pends on both the graph of the system and the strengths of the 
interactions in the Hamiltonian. 

We speculate that the revival structure described in the 
previous sections is much more generic and valid whenever 
a many-body quantum system has local dynamics, as long 
as the integrability breaking is not too strong and hence 
information is not completely lost in the system. Since 
the propagation of information is governed by vlr, the 
revival time scale becomes T rov ~ L/vlr in general. In the 
integrable case, vlr reduces to the quasiparticle speed v max , 
and hence we recover the revival time scale in Eq. ([8]). On 
the other hand, the locality of the dynamics is intact in the 
non-integrable case as well, providing a natural explanation 
for the robustness of the revival structure. 



VI. SUMMARY 

In this paper, we studied the phenomenon of revivals after a 
quantum quench in the transverse field XY model, and found 
a non-trivial revival structure that can not be obtained from 
a simple spectral analysis. It was shown that this structure 
is universal in the sense that it does not depend on the initial 
state and the size of the quench. Revivals were shown to be 
related to quasiparticles propagating around the system with a 
finite maximum speed i> max , and a corresponding estimate for 
the revival time scale was established. 

We also investigated the effect of non-integrable perturba- 
tions on the structure of the revivals. In particular, we consid- 
ered the XZ spin chain and a random disorder in the magnetic 
field. It was found that the revival structure is clearly visible at 
time scales far beyond the lifetime of quasiparticles, implying 
that something more generic than integrability is behind the 
phenomenon of revivals. In perspective of this, we proposed 
a generic connection between revivals and locality, where the 
quasiparticle speed v max generalizes into the Lieb— Robinson 
speed Vlr- We believe that a thorough understanding of this 
important connection requires further study of non-integrable 
systems, for example, studying how the entanglement produc- 
tion in the subsystem is related to the loss of the revival struc- 
ture. 
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